Many-Body Electrostatic Forces Between Colloidal Particles at Vanishing Ionic 
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Electrostatic forces between small groups of colloidal particles are measured using blinking optical 
tweezers. When the electrostatic screening length is longer than the inter-particle separation, forces 
are found to be non-pairwise additive. Both pair and multi-particle forces are well described by the 
linearized Poisson-Boltzmann equation with constant potential boundary conditions. These findings 
may play an important role in understanding the structure and stability of a wide variety of systems, 
from micron-sized particles in oil to aqueous nanocolloids. 



Colloidal suspensions offer a unique opportunity to 
observe the relationship between microscopic interac- 
tions and macroscopic thermodynamic behavior because 
the component particles are both large enough to ob- 
serve directly with light microscopy and small enough to 
come to thermodynamic equilibrium over experimental 
timescales. Colloids display a rich range of phase behav- 
ior including equilibrium gases, liquids, and crystals and 
nonequilibrium glasses and gels. Controlling and predict- 
ing this phase behavior has applications from ensuring 
the stability of paint to directing the self assembly of 
novel optical materials. 

Predicting phase behavior requires an understanding 
of the interactions of the component particles. Forces 
between pairs of charged surfaces have been extensively 
characterized using direct force measurements Al- 
ternatively, observations of bulk correlations can be used 
to extract an "effective pair potential," but this pair po- 
tential is found to vary with the density of the suspen- 
sion 0, 0| ■ This suggests that charged colloidal particles 
do not, in fact, interact only through a pair potential. 
Many-body interactions evidently play an important role 
in determining macroscopic behavior. 

To bridge the gap between pair interactions and bulk 
behavior, we explore the regime of few-body interac- 
tions by using blinking optical tweezers to directly mea- 
sure electrostatic forces between sets of several isolated 
charged colloidal particles at low ionic strength. For 
screening lengths much shorter than the inter-particle 
separation, the measured forces are nearly pairwise addi- 
tive; however, for screening lengths longer than the inter- 
particle separation, there are significant deviations from 
pairwise additivity. 

Our system consists of radius a = 600 nm PMMA 
spheres suspended in hexadecane, a non-polar solvent. 
The spheres are coated in PHSA in order to stabilize 
them against aggregation Both the particle charge 
and screening length are adjusted by adding variable 
concentrations of a reverse micelle forming surfactant, 
NaAOT, as described elsewhere 0,(1,0,11. 



Forces are extracted from the trajectories of freely in- 
teracting particles, which are repeatedly trapped and re- 
leased using blinking optical tweezers [9( created with a 
1064 nm laser. The particles are more than 30 particle 
diameters from the nearest wall. A detailed explanation 
of the apparatus is given in [loj | . The physical ideas un- 
derlying the technique are developed in 11[ . The central 
result of that paper is a relation between the force, /, and 
the two statistical parameters that describe ensembles of 
short-time Brownian trajectories, the drift velocity, Vd, 
and the diffusion coefficient, D: 



f = ksTD Vd- 



(1) 



This relation holds for a single particle in one dimension, 
or separately for each of the hydrodynamic normal modes 
of a multi-particle system in several dimensions. The 



range of validity of this expression is explored in 11|. 

We measure forces between beads in pair, equilateral 
triangle, and hexagonal configurations, as depicted in 
Fig. [TJ The first column of Fig. [TJ shows the measured 
force as a function of separation, s, between pairs of beads 
at two different surfactant concentrations: [AOT] = 10 
mM (top) and 0.5 mM (bottom). At each separation, 
three beads are chosen for pair measurements. The forces 
between each of the three possible pairs formed from this 
set are plotted with different colors. For a given concen- 
tration and separation, forces between each of the pairs 
are nearly identical, which suggests that each of the beads 
has nearly the same charge. The force at the smallest 
measured separation is roughly the same for both con- 
centrations, but the interaction is screened much more 
rapidly at the higher surfactant concentration. 

Next, the same three beads are formed into an equi- 
lateral triangle. The forces on the breathing mode are 
calculated using Eq. [TJ and displayed in the second col- 
umn of Fig. [TJ The previously measured pair forces are 
used to make a pairwise prediction of the force on the 
equilateral triangle. Since the same beads are measured 
at the same separations for both the pair and equilat- 
eral geometry, this experiment provides a strong, model- 
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FIG. 1: Direct Measurement of Non-Pairwise Electrostatic Interactions. Forces on beads in pair (first column), equilateral 
(second column), and hexagonal (third column) configurations, at AOT concentrations of 10 mM (first row) and 0.5 mM 
(second row). The arrows on the particle configurations in the first row indicate the form of the breathing modes used for 
analysis. Breathing modes are normalized so that sum of the squares of the particle displacements is unity. For the pair 
measurements, different colored points represent different pairs in the same sample. The dotted, solid, and dashed lines are 
fits to constant charge density, constant potential, and a simple approximation of constant potential based on Eq. [3] For the 
equilateral and hexagonal configurations, black points are measured forces on the breathing mode and red points are a direct 
pairwise sum of the measured pair forces. The red lines are pairwise sums of the constant potential pair fits. Constant potential 
predictions for the force on the breathing mode based on fits to the pair data are shown as black lines. The solid line is based 
on the full numerical solution while the dashed line is based on Eq. [3] 



independent test of pairwise additivity. At 10 mM, agree- 
ment with the pairwise prediction is very good, whereas 
at 0.5 mM, there are significant deviations from pairwise 
additivity. 

Finally, the same three beads along with four others 
are formed into a hexagon with a particle at its center. 
The forces on the breathing mode are displayed in the 
third column of Fig. [1] There is a striking disagreement 
with the pairwise additive prediction at 0.5 mM, and a 
small deviation at 10 mM. 

At the end of the experiment, after measurements in all 
three configurations are complete, force measurements at 
the smallest separation in the pair and hexagonal config- 
urations are repeated. These measurements are plotted 
as open symbols. Their agreement with measurements 
from the beginning of the experiment suggests that the 
forces do not drift noticeably over time. 

These forces are electrostatic. In the presence of mobile 
ions, the electrostatic potential, ip, obeys the Poisson- 



Boltzmann equation (PBE): 

/3eV 2 G- = K 2 smh(/3eV0, ( 2 ) 

where k 2 — 2n oc f3e 2 / e is the ionic strength, which is pro- 
portional to the density of ions in the bulk, n.^. The 
quantity k -1 , called the screening length, sets the range 
of electrostatic interactions. Since the screening length 
scales , the dimensionless parameter na 

can be the same for micron-sized particles in oil and 
nanocolloids in water. 

As a result of the non-linearity of the PBE, its solu- 
tions cannot be superposed. Thus, forces are not gener- 
ally pairwise additive. Russ et al. numerically solved the 
PBE with constant charge density boundary conditions 
for two and three spheres to find a three-body potential 
[ 1 j . This theory was used to explain deviations from 
pairwise additivity found in experiments conducted in 
water with kcl « 2 [13}. The constant charge density 
PBE fails to describe our data, where na < 1, because it 
predicts deviations from pairwise additivity that increase 
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with ionic strength: the opposite of the trend observed 
in our system. 
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FIG. 2: Charge reduction in the presence of a third particle 
Ratio of the charge per particle on a set of three constant 
potential spheres arranged in an equilateral triangle of side 
length s to the charge per particle on a pair of similar spheres 
separated by length s for six values of the ionic strength: na — 
0, 0.2, 0.4, 0.6, 0.8 and 1.0. The solid line represents a nu- 
merical solution of the LPBE, and the dashed line represents 
an analytical approximation based on a single point charge at 
the center of each sphere, Eqs.[UE] 

Non-linear governing equations are not necessary for 
non-pairwise additive forces. Indeed, forces between fi- 
nite size conductors held at constant potential are non- 
pairwise additive even in free-space electrostatics, which 
is a perfectly linear theory. In fact, we find that the lin- 
earized Poisson-Boltzmann equation (LPBE) is sufficient 
for describing all of our observations including deviations 
from pairwise additivity. 

In order to maintain constant potential, the total 
amount of charge on a pair of spheres must change as 
they are brought together. As a result, the force law 
between a pair of conducting spheres is different from 
the force law between point charges 



14j . Similarly, as 



a third sphere approaches a pair, the charge will change 
even further. This is the reason for non-pairwise additiv- 
ity of forces. This is illustrated in Fig. [2j where we plot 
the ratio of the charge on each sphere in an equilateral 
triangle of side length s to the charge on each sphere in 
a pair of the same separation. The charge is calculated 
for several values of na by numerically solving the LPBE 
using a simple singularity method [15j ]. 

Fits to the constant potential pair force law 
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shown as solid lines in Fig. [TJ The resulting values of na 
and ?po are shown in Fig. [TJ For comparison, fits to the 
constant charge density force law are shown as dotted 
lines. Clearly, it is impossible to distinguish between the 
two force laws based only on these pair force data. How- 
ever, the two fits give different values for the screening 
length and the charge or surface potential. This leaves 
an essential ambiguity about the underlying physics. 
The many body data resolve this ambiguity. Based 



on the surface potential and screening length extracted 
from the constant potential fits, we calculate the ex- 
pected forces on the breathing mode of the equilateral 
and hexagonal configurations, which are plotted as solid 
black lines in Fig.[TJ This linear constant potential model 
gives excellent quantitative predictions for the deviations 
from pairwise additivity observed at [AOT] = 0.5 mM. 
The constant charge density model predicts no devia- 
tions from pairwise additivity. At 10 mM, where both 
predicted and observed deviations from pairwise additiv- 
ity are small, it is not clear that the constant potential 
model is much better than simply assuming pairwise ad- 
ditivity. We emphasize that all the information used to 
calculate the expected forces in the equilateral and hexag- 
onal configurations comes from fits to the pair data — no 
additional free parameters arc introduced. 

In order to understand why this constant potential 
model predicts larger deviations from pairwise additivity 
at lower ionic strength, it is helpful to consider a sim- 
ple analytical approximation. Model each particle with a 
single point charge at its center with a value that is free 
to change in order to make the particle surfaces roughly 
constant potential. Solving the LPBE shows that the po- 
tential on the surface of a sphere due to a point charge q 
at its center is qe~ Ka /Airea. The potential due to a point 
charge outside the sphere, a distance s from its center, 
varies over the surface of the sphere, but can be roughly 
approximated by the value of the potential at its center, 
qe~ KS /Aires, so long as s > o. The total potential on the 
surface of particle i is thus approximately 



V'o.i = -7— ( H 6 
Aire \ a 



(3) 



where is the separation between particles i and j. If 
we demand that this potential takes on some fixed value, 
then this is a system of linear equations that can easily 
be solved for the charge at the center of each particle as 
a function of the configuration, s,j . The solutions for the 
charge per particle in pair and equilateral configurations 
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AneM^^/a + e^/s)- 1 (4) 
4neiP (e~ Ka /a + 2e- KS /s)- 1 (5) 



It is important to note that the singular charge, q, is 
only equal to the actual charge on a particle, Q, in the 
Coulomb limit. In general, for isolated particles they 
are related by Q = qe~ Ka (l + na), a result that can be 
obtained from Gauss's law. Predictions based on this 
approximation for the total charge as a function of sepa- 
ration are plotted as dashed lines in Fig. O The approx- 
imation is apparently quite good. 

In the equilateral configuration, the force on any one 
of the spheres should scale as q 2 , so the ratio of the true 
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force to the pairwise force is, in this approximation, 
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identical constant potential particles as a function of sep- 

At very low 
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(6) 

This shows how the non-pairwise additivity depends on 
the three relevant length scales: s, and a. The ex- 

pected forces based on this approximation are plotted 
as dashed lines in the second column of Fig. [1] and are 
in excellent agreement with both the numerical calcula- 
tion and the data over the full range of measured sep- 
arations. This analytical approximation does not take 
polarization into account, whereas the numerical calcu- 
lation does. Their agreement demonstrates that, in these 
experiments, polarization is not important compared to 
the change in total charge. 

Assuming the surface potential is independent of con- 
figuration leads to an accurate prediction of the observed 
deviations from pairwise additivity. Assuming constant 
charge density is evidently wrong, at least for the mea- 
surements at lower ionic strength. But why should the 
surface potential be constant? After all, these are iso- 
lated dielectric particles, not conductors hooked up to a 
battery. 

In fact, the idea of constant potential colloid surfaces 
is far from new. In their foundational 1948 work on col- 
loidal interactions, Verwey and Overbeek argued that 
colloidal particles should be at constant surface poten- 
tial whenever surface ions are in equilibrium with bulk 
ions, and criticized earlier assumptions of constant charge 
as artificial [lij]. Finding equilibrium by equating the 
bulk and surface ion chemical potentials led them to the 
Nernst equation, /3c0o = hi(rioo/«-rcf ), where n re f is a 
material parameter. This suggests that the surface po- 
tential depends only on material properties and the bulk 
concentration of ions, and not on geometry or any inter- 
actions between particles. However, they explicitly ig- 
nore the surface configurational entropy, which can cou- 
ple the surface potential to the fields to push the system 
toward constant charge 17, 18, Different models of 
the surface chemistry make different predictions for the 
surface configurational entropy, so the agreement of con- 
stant potential boundary conditions with our data places 
constraints on the surface chemistry. In particular, it ap- 
pears that a surface with both positive and negative dis- 
sociable ions — an amphoteric surface — should be in bet- 
ter accord with the Nernst equation and thus closer to 
constant potential than a surface with only a single dis- 
sociable ion species [20I . 



With few exceptions [21| . previous predictions of col- 
loidal structure have uniformly assumed constant charge 
density boundary conditions @, [SJ 



23, |2J, |25(. Our 



experimental evidence for constant potential boundary 
conditions demands a fresh look at the bulk structure of 
charged colloidal suspensions. As a first step, we have 
calculated the charge per particle in a bulk system of 



aration using the PB cell model 2d 
ionic strength, the charge per particle drops significantly 
even at modest volume fractions. For na — 0.14 and 
(ie^o = 3.4, the measured values at [AOT] = 0.5 mM, 
the charge per particle at 10% volume fraction is reduced 
by a factor of 6 compared to an isolated particle. Thus, 
we expect constant potential boundary conditions to have 
important implications for predictions of aggregation and 
self assembly in a wide range of systems where the screen- 
ing length is larger than the particle size. An interesting 
case to consider is binary oppositely charged colloids at 
kcl « 1, where the self assembly of diamond- like struc- 
tures has been observed [28J. In that case, constant po- 
tential boundary conditions should have the opposite ef- 
fect on particle charge — as the volume fraction increases 
we expect the charge per particle to increase. Much the- 
oretical work is needed. 
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